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Abstract. The study of the B + — > ppK + decay channel at LHCb offers great opportunities to study different 
aspects of the Standard Model and possibly Beyond Standard Model physics. In particular it can be interesting not 
only for the possibility to measure CP asymmetry but also to study possible intermediate resonances. The ratios of 
the branching fractions of the B + — > ppK + decay channel, of the charmless component with M p p < 2.85 GeV/c 2 
and of the charmonium contribution rj c relative to the J/tf/ are presented. 



1 Introduction 

The current knowledge on the B + — > ppK + decay is based 
on the measurements performed at the /^-factories II 1 121 . 
Thanks to the high statistics available in the next years 
at LHCb, due to the large bb production cross section, 
cr h - b = (284 + 20 + 49)/i b Q at a centre of mass energy 
of -\fs = 7 TeV, more precise measurements can be per- 
formed and possibly new physics can be revealed. 
LHCb is a single arm forward spectrometer optimised to 
perform precise CP violation measurements and rare de- 
cays study in the heavy flavour sector. 
Three-body baryonic B decays are expected to proceed 
predominantly via b — > s penguin diagrams or doubly 
Cabibbo suppressed tree diagrams |4| and as such new 
heavy virtual particles can give measurable departures 
from Standard Model expectations. The baryon-antibaryon 
pair can also be produced by a pair of gluons in OZI 
suppressed penguin diagrams, where gluonic resonances 
could be formed. The three Feynman diagrams are shown 
in Fig. [T] The interference between the penguin diagram 




Fig. 1. Feynman Diagrams: penguin (top left), doubly Cabibbo 
suppressed tree (top right), OZI suppressed penguin (bottom). 



and the doubly CKM suppressed tree diagram can lead to 
a direct CP asymmetry. 
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Moreover, three-body baryonic decays offer a clean en- 
vironment to study intermediate states, as charmonium 
states, excited A baryons and exotic states (glueballs Q, 
baryonium, pentaquarks). In fact, the pp final state can be 
produced through a resonant intermediate state, a cc state 
or through an unknown exotic state. At LHCb, beyond the 
study of well known charmonium states, it can be interest- 
ing to observe charmonium-like states such as X(3872) or 
new intermediate states in order to study their nature and 
properties. 

Another interesting feature observed both at the BaBar and 
Belle detectors in three body baryonic B decays is the low 
mass M p p enhancement. Several explanations have been 
proposed. A study of the distribution of events in the Dalitz 
plot can explain the origin of this mass enhancement. j6). 
A preliminary study of this mode based on the sample 
of events collected by the LHCb detector in the 2010 
data taking period is presented. The branching ratios of 
B + — > ppK + , of the charmless component in the range 
Mpp < 2.85 GeV/c 2 and of the B ± -> ricK* decay channel 
with r\ c — > pp relative to the product of branching fractions 
S m = SiB* -> J/1//K*) x S(J/t// -> pp), which is well 
known [7], have been measured. 

The measurement relies only on the ratios of events, 
Af(mode), and efficiencies, e mo d e , with respect to the ref- 
erence mode: 

S(mode) JV(mode) e, /lA 

— X (1) 

where N(J/i/r) = N B+ ^j/^p p)K+ ande y/(A = e s ^;/^ P p )r 

2 Study of the B + -> ppK + selection 

Since the bb cross section is approximately 1 % of the total 
inelastic cross section and the interesting channel branch- 
ing ratios are of the order of S ~ 10 5 -10 6 , it is extremely 
important to design a selection with a high rejection of the 
background. The selection strategy is based on a preselec- 
tion with very loose selection criteria. A multivariate anal- 
ysis using the "Boosted Decision Tree" algorithm of the 
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Fig. 2. Mass distribution of all selected B + — » ppK + events. The 
continuous line shows the fit result, the dashed line shows the 
fitted background component. 
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Fig. 3. Mass distribution of the selected fi + — > ppK + events for 
M p p < 2.85 GeV/c 2 . The continuous line shows the fit result, the 
dashed line shows the fitted background component. 



TMVA package [8| is used to optimise the selection effi- 
ciency. The variables used in the selection procedure ex- 
ploit kinematic and topological characteristics of fe-hadron 
decays. Three charged particles in the final state with large 
transverse momenta are required. The three final state par- 
ticles have to come from the same displaced secondary ver- 
tex requiring a good vertex fit and a significant distance 
between the two vertices. Moreover it is required a small 
maximum distance of closest approach between the tracks 
and non zero impact parameters. And the momenta direc- 
tion of the three daughters has to be compatible with the 
flight direction of the B. Furthermore the particle identi- 
fication information, provided by two RICH detectors in 
LHCb, is a fundamental ingredient of the selection in or- 
der to identify correctly protons and kaons candidates. A 
signal efficiency of more than 70% with respect to prese- 
lected events has been obtained using the BDT algorithm 
with a signal to background ratio of S /B ~ 1. 



3 Signal Extraction 

The number of B + — > ppK + signal events is extracted from 
an unbinned maximum likelihood fit to the M ppK invari- 
ant mass distribution of the selected events (Fig. El. The 
fit procedure has been applied also on the subsample of 
events having M pp < 2.85 GeV/c 2 to extract the number of 
events in the charmless region (Fig. [3]). The signal proba- 
bility density function has been parametrized with a Gaus- 
sian function while for the background a linear function 
has been used. 

For the B + -> J/if/(-^> pp)K + and B + /? e (-> pp)K + , the 
yield has been extracted with a fit to the M pp distribution 
for the B events having M pp k invariant mass compatible 
within 2.5<t with the known B + mass (Fig. |4j). Since the 
J/if/ meson has a natural width much smaller than the de- 
tector resolution while the natural width of the t] c is not 
negligible, the J/t// has been parametrized with a Gaussian 
function while for the r\ c the convolution of a Breit-Wigner 
with a Gaussian function with the same resolution of the 
J/if/ has been used. 



4 Estimation of efficiencies 

The calculation of the branching fractions based on Eq. [T] 
requires the determination of the ratio of the efficiency 
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Fig. 4. Mass distribution of pp in the J/if/-r] c range for events 
having \M(ppK) - M(B)\ < 50MeV/c 2 . The solid line shows the 
result of the fit, the dashed line corresponds to the fitted back- 
ground contribution. 

for each considered channel (all the events, events with 
M pp < 2.85GeV/c 2 and j] c events) with respect to the 

The overall efficiency ratios, product of reconstruction, trig- 
ger and event selection efficiency as a function of M pp in- 
variant mass have been calculated using simulated events 
in order to correct for varying efficiency over the M pp range. 
The ratios with their uncertainties are reported in Tab. [T] 

5 Evaluation of systematics 

Since the final state for the considered modes and the refer- 
ence one is the same, most of the systematic uncertainties 
cancel. 

Possible systematic uncertainties have been considered on 
the event yield determination and on the efficiency ratios. 
The systematic uncertainty on the event yield extraction is 
estimated using different fit ranges, different parametrisa- 
tion for signal and background. Half of the largest differ- 
ence in the central value obtained with different fit config- 
urations is taken as systematic uncertainty on the number 
of signal events. The event yields for each mode and the 
corresponding systematic uncertainties are listed in Tab. [T] 
The uncertainty on the efficiency ratios described in Sec- 
tionals considered as a systematic uncertainty. 
Most data-simulation discrepancies are independent from 
the Dalitz plot variables and cancel in the ratio. In some 
cases, in particular for particle identification, the discrep- 
ancies are larger for low momentum tracks. A conserva- 
tive 10% systematic uncertainty due to momentum depen- 
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dent data-simulation discrepancy on the efficiency ratios 
has been assumed for the total T>(B + — > ppK + ) and for the 
charmless S(B + — > ppK + ) Mpj . <2 ^Q e y/ c i branching ratios. 



6 Branching ratio measurements 

Using the number of events extracted for each considered 
mode and the corresponding ratios of efficiencies with re- 
spect to the B + — > 7/i/'(—> pp)K + the preliminary measure- 
ments of ratios of branching fractions have been extracted 
and are reported in Tab. [2] 

Table 1. Yields of B + — > ppK + events in each mode, ratio of ef- 
ficiencies of each mode with respect to the B + — » J/i//{— » pp)K + . 



Mode 


Yield 






± stat ± syst 






48 ± 7 ± 2 


1 


ALL 


225 ± 18 ± 5 


1.01 ±0.10 


Mpp < 2.85 GeV/c 2 


102 ± 12 ±2 


0.96 ±0.10 


1c 


34 ± 8 ± 3 


0.994 ± 0.003 



Table 2. Ratios of branching fractions with statistical and sys- 
tematic uncertainties for the different modes considered. 



Mode 


S ratio 




LHCb 




1 


ALL 


4.6 ±0.6 ±0.5 


M pr , < 2.85 GeV/c 2 


2.21 ±0.41 ±0.24 


1c 


0.71 ±0.20 ±0.07 



7 Conclusions and perspectives with the 
full 201 1 data set 

The preliminary measurements of the ratios of branching 
fractions for the B + — > ppK + decays have been reported. 
From the measured ratios, using the known values of the 
S(B + -> J/ipK + ) and S(J/i// -> pp) the corresponding 
branching fractions have been derived and reported in Tab. [3] 



Table 4. Branching ratios measured at the S-factories. All values 
are in units 10~ 6 . 



Mode 


BaBar 


Belle 




m 


□ 


ALL 




10.76!^ ± 0.70 


Mp P < 2.85 GeV/c 2 


5.3 ± 0.4 ± 0.3 


5.0^ ± 0.32 


»7c 


1.8^ ±0.2 


1.42 ± 0.11^ 



The results are compatible with the average values re- 
ported in the PDG and measured at the fi-factories as re- 
ported in Tab. [4] even with the limited statistics used in this 
analysis. 

The current systematic uncertainty is largely of statis- 
tical origin, so both statistic and systematic uncertainties 
will decrease as soon as the available statistics increase. 

In 201 1 LHCb has collected more than 1 ftr 1 . The large 
available statistics collected in 2011 is a factor almost 10 
more with respect to the statistics collected at the fi-factories 
The most accurate available measurements have been per- 
formed by the BaBar and Belle Collaborations using about 
1000 fully reconstructed exclusive decays. An improve- 
ment on the precision of the branching fractions is ex- 
pected in particular for < B(B ± — > rjcK 12 ) x £(t] c — > pp). 
Moreover the study of more rare charmonium contribu- 
tions will be possible, such as r] c (2S) and X(3872), observ- 
ing or putting an upper limit on S(B + — > ccK + ) x < B(cc — > 

pp) ~ io- 7 . 
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Table 3. Derived LHCb preliminary branching ratios with sta- 
tistical error, systematic uncertainty and the uncertainty on the 
knowledge of the S(B + — > J/ifr(—> pp)K + ). All values are in 
units 10~ s . 



Mode 


LHCb S 
(preliminary) 


ALL 


10.2 ± 1.4 ± 1.1 ±0.5 


Mp r , < 2.85 GeV/c 2 


4.87 ±0.91 ±0.54 ±0.22 


1c 


1.57 ±0.43 ±0.15 ±0.07 



